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ABSTRACT 

Recent quantum mechanical calculations of rate coefficients for collisional transfer 
of population between the hyperfine states of enable their population densities 

to be determined. We have computed the relative populations of the hyperfine states 
of the N = 0,1,2 rotational states for kinetic temperatures 5 ^ T ^ 20 K and 
molecular hydrogen densities 1 ^ n(H 2 ) ^ 10^° cm“^. Spontaneous and induced 
radiative transitions were taken into account. Our calculations show that, if the lines 
are optically thin, the populations of the hyperfine states, F, within a given rotational 
manifold are proportional to their statistical weights, {2F + 1) - i.e. in LTE - over the 
entire range of densities. We have re-analyzed IRAM 30 m telescope observations of 
hyperfine transitions {N = 1 —>■ 0) in four starless cores. A comparison of these 
observations with our calculations confirms that the hyperfine states are statistically 
populated in these sources. 
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1 INTRODUCTION 

The interpretation of the emission lines of molecules in the 
interstellar medium (ISM) is often complicated by the ef¬ 
fects of re-absorption and scattering, owing to significant 
optical depths in the lines. Partly for this reason, observa¬ 
tions of less abundant isotopologues are analyzed, in addi¬ 
tion or in preference to those of the principal species; this 
is the case of CO, for example, where and C^®0 lines 

are used, and also of CN, where and serve a 

similar purpose. In the present paper, we consider the emis¬ 
sion lines of ^®CN, observed at millimetre wavelengths in 
pre-protostellar sources. 

has a rich spectrum at mm-wavelengths, where 
it displays the effects of the fine structure interaction, be¬ 
tween the electron spin and the nuclear rotation, and of the 
hyperfine interaction with the spins of the {h = |) and 
N (U = 1) nuclei. Thus, Bogey et al. (1984) listed 16 ‘al¬ 
lowed’ transitions, in the vicinity of 100 GHz, between the 
rotational states N = 0 and N = 1, and a further 25 transi¬ 
tions at 200 GHz between the N = 1 and N = 2 states. It is 
this cornucopia of optically thin transitions that one wishes 
to exploit, in order to obtain a better understanding of the 
conditions in pre-protostellar objects. 

It is generally assumed that the hyperfine levels, F, of 
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a given N are populated in proportion to their statistical 
weights, {2F -\- 1), i.e. that they are in local thermodynamic 
equilibrium LTE. However, LTE is the exception, rather 
than the rule, in the ISM, because very low densities pre¬ 
vail. The assumption of LTE is usually dictated by a lack of 
rate coefficients for collisional population transfer between 
the hyperfine levels of the molecule. In the case of ^®GN 
(and of C^®N), this situation has been rectified recently by 
the calculations of Flower & Lique (2015), which provide 
rate coefficients for collisions with para-H 2 in its rotational 
ground state - the dominant perturber at low kinetic tem¬ 
peratures, T. Thus, the opportunity arises to calculate the 
hyperfine level populations explicitly, allowing for collisional 
and radiative transfer, where the latter is induced by the 
cosmic microwave background or by the emission of dust 
present in the medium. 

We present the details of the calculations in Section 2 
and the numerical results and their implications in the fol¬ 
lowing Section 3. Our concluding remarks are in Section 4. 


2 CALCULATIONS 

Bogey et al. (1984) note that the rotational (A) and the 
hyperfine {F) are the ‘good’ quantum numbers of ^^CN. 
Accordingly, we evaluate the relative populations, in equi¬ 
librium, of the levels {N,F), for a wide range of values of 
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the density of molecular hydrogen (specifically, para-H 2 ), 
1 ^ n(H 2 ) ^ 10^° cm“®, and a range of kinetic tempera¬ 
tures, 5 ^ T ^ 20 K. Collisional and radiative population 
transfer are taken into account, as described below. 


2.1 Collisional population transfer 

We use the rate coefficients for collisions between and 
para-H 2 that were calculated by Flower & Lique (2015), who 
used the angular momentum coupling scheme j = N + S, I = 
Ii +12, F = j + 1 , where N is the rotational quantum num¬ 
ber, S' = i is the net electron spin, Ii = \ and /2 = 1 are the 
net spin quantum numbers of the and N nuclei, respec¬ 
tively, and F is the hyperfine quantum number. Their calcu¬ 
lations were performed separately for the cases of 7 = ^ and 
/ = |, as the value of the total nuclear spin, 7, is conserved 
during the collision. For a given value of 7, the rate (s“^) 
of collisional transfer of population from the state {N, j, F) 
to the state {N', j', F') is n{H 2 )<}{N, j, F —>■ N', j', F'), and 
hence the total rate of transfer from level (N, F) to level 
{N',F') is n(H 2 )^,_,,g(iV,j,F ^ N',j',F'), where q de¬ 
notes a rate coefficient (cm^ s“^) and n(H 2 ) is the density 
(cm“®) of ground state para-H 2 . We assume that the nuclear 
spin states 7 = i and 7 = | are statistically populated and 
evaluate the rate coefficients as 

q{N,j,F^N',j',F') = ig(7= i|7V,i,F^ 7V',j',F')+ 

lq{I^l\N,j,F^N',j',F'). (1) 

Varying the relative populations of the nuclear spin states 
had no significant effect on the results of the calculations 
reported in Section 3 below. 


2.2 Radiative population transfer 

Following the analysis of optically allowed transitions by 
Nicholls & Stewart (1962), the rate (s“^) of spontaneous 
radiative transitions between states i and j, belonging to 
levels k and I, respectively, with energies Ek < Ei, is given 
by 


(V, Fi , F 2 , F), where Fi = S-f-Ii, F 2 = N-I-Fi, F = F 2 -fl 2 
is the angular momentum coupling scheme adopted by Bo¬ 
gey et al. (1984). The A-values for transitions between pairs 
of these states are given in Table 1. The accuracy of these 
results is determined principally by the accuracy of the 
dipole moment, which is difficult to assess from the available 
data (http://cccbdb.nist.gov/dipole2.asp). The rates (s~^) 
of spontaneous transitions between levels {N, F) were ob¬ 
tained by implementing the summation relation (4) above. 

According to the Einstein relations, the rates of radia¬ 
tive excitation and de-excitation by a radiation field of en¬ 
ergy density p^du (erg cm“®) at frequency u are Bkipv and 
Bikpv, respectively, where 


^kBki = ijJiBik = 


STThl/^ 


uiAik 


(5) 


and bj = (2F-I-1) is the statistical weight. We adopt a black- 
body radiation field. 


Pi/( 7b) 


1 

'''' exp - 1 


(6) 


where Th = 2.73 K is the temperature of the cosmic back¬ 
ground radiation (CBR) or of a black-body radiation field 
produced by local dust, mixed with the gas. 


2.3 Level populations and line intensities 

Knowing the total rates. 


Pki = n{R2)qki{T) + Aki + Bki, (7) 

of collisional and radiative population transfer between lev¬ 
els {N,F), the populations, rik, of the 12 levels with rota¬ 
tional quantum numbers N = 0,1, 2 may be obtained by 
solving the rate equations 

= y^^jmpik - rikPki) ( 8 ) 

l^k 

together with the closure relation 


Aji — 


4(Fi - Ekf 


where the line strength 


( 2 ) 


Sji — ^jip 5 (3) 

and where p = 1.45 debye (cf. Hily-Blant et al. 2013) is 
the electric dipole moment of the molecule in its vibrational 
ground state. The algebraic quantities sji were computed by 
Bogey et al. (1984, table 3) for dipole transitions within the 
V = 0 ■4— 1 and V = 1 4— 2 rotational manifolds. We note 
that 


Y,nk=nC^CN), (9) 

k 

where n(^^CN) is the number density of molecules. In 

steady-state, = 0, and the relative level populations are 
readily obtained by matrix inversion, setting n(^®CN) = 1. 
Then, the relative line emission (erg s“^) is given by 

Iji = niAjihuji, ( 10 ) 

where, we recall, j and i are states {N, F\, F 2 , F) belonging 
to the levels I and k, identified by the quantum numbers 
{N,F); 12 ji is the line frequency. 


ij 2.4 Line opacities 


In the present context, the levels k, I are identified by The opacity at frequency u in a transition of oscillator 

the quantum numbers {N,F), and the states i,j by strength/ij is given by 
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Table 1. Spontaneous radiative transitions probabilities, A, and relative line emission, I [equ. (10)], for transitions = 0 N' = 1 and 
N = 1 <— N' = 2; T = 10 K and n(H 2 ) = 10^ cm“^. The angular momentum coupling scheme of Bogey et al. (1984) is adopted. The 
quantity s is the algebraic factor in the expression for the line strength, equ. (3), and u is the line frequency. Note that the degeneracy 
of the upper (emitting) level is {2F' + 1). Numbers in parentheses are powers of 10. 


N 

Fi 

F2^F' 

F <- F' 

s 

u MHz 

A s-i 

I erg s ^ 

0 

0 

0 


1 

1 


2 

1.66 

108651.297 

1.04(-5) 

1.39(-1) 

0 

0 

0 


1 

1 


1 

1.00 

108636.923 

1.05(-5) 

8.57^2) 

0 

0 

0 


1 

1 


0 

0.33 

108631.121 

1.04(-5) 

2.81(-2) 

0 

1 

1 


2 

0 


1 

0.56 

108786.982 

5.88(-6) 

4.80(-2) 

0 

1 

1 


2 

1 


2 

1.25 

108782.374 

7.87(-6) 

1.05(-1) 

0 

1 

1 


2 

2 


3 

2.33 

108780.201 

1.05(-5) 

1.94(-1) 

0 

1 

1 


2 

1 


1 

0.42 

108793.753 

4.41(-6) 

3.60(-2) 

0 

1 

1 


2 

2 


2 

0.42 

108796.400 

2.65(-6) 

3.54(-2) 

0 

1 

1 


1 

0 


1 

0.33 

108638.212 

3.45(-6) 

2.81(-2) 

0 

1 

1 


1 

1 


2 

0.42 

108643.590 

2.64(-6) 

3.53^2) 

0 

1 

1 


1 

1 


1 

0.25 

108645.064 

2.62(-6) 

2.14(-2) 

0 

1 

1 


1 

2 


2 

1.25 

108657.646 

7.85(-6) 

1.05(-1) 

0 

1 

1 

<— 

1 

1 


0 

0.33 

108644.346 

1.04(-5) 

2.81(-2) 

0 

1 

1 


1 

2 


1 

0.42 

108658.948 

4.39(-6) 

3.58^2) 

0 

1 

1 


0 

1 


1 

0.33 

108412.862 

3.43(-6) 

2.79(-2) 

0 

1 

1 


0 

2 

4 - 

1 

0.56 

108426.889 

5.82(-6) 

4.74(-2) 

1 

0 

1 


2 

0 


1 

0.33 

217296.605 

2.76(-5) 

3.00(-2) 

1 

0 

1 


2 

1 


2 

0.75 

217301.175 

3.77(-5) 

6.25(-2) 

1 

0 

1 


2 

2 


3 

1.40 

217303.191 

5.02(-5) 

1.19(-1) 

1 

0 

1 


2 

1 


1 

0.25 

217290.823 

2.09(-5) 

2.27(-2) 

1 

0 

1 


2 

2 


2 

0.25 

217286.804 

1.26(-5) 

2.09(-2) 

1 

1 

2 


3 

1 


2 

0.84 

217469.151 

4.23(-5) 

7.02(-2) 

1 

1 

2 


3 

2 


3 

1.24 

217467.150 

4.46(-5) 

1.06^1) 

1 

1 

2 


3 

3 


4 

1.80 

217467.150 

5.03(-5) 

1.61(-1) 

1 

1 

2 


3 

2 


2 

0.16 

217480.559 

8.05(-6) 

1.34(-2) 

1 

1 

2 


3 

3 


3 

0.16 

217483.606 

5.75(-6) 

1.37(-2) 

1 

1 

1 


2 

0 


1 

0.25 

217443.722 

2.10(-5) 

2.29^2) 

1 

1 

1 


2 

1 


2 

0.57 

217436.350 

2.87(-5) 

4.76(-2) 

1 

1 

1 


2 

2 


3 

1.05 

217428.563 

3.77(-5) 

8.96^2) 

1 

1 

1 


2 

1 


1 

0.19 

217443.722 

1.59(-5) 

1.73(-2) 

1 

1 

1 


2 

2 


2 

0.19 

217437.702 

9.56(-6) 

1.59^2) 

1 

1 

2 


2 

1 


1 

0.11 

217294.470 

9.21(-6) 

1.00(-2) 

1 

1 

2 


2 

2 


2 

0.17 

217298.937 

8.54(-6) 

1.42(-2) 

1 

1 

2 


2 

3 


3 

0.31 

217306.117 

l.ll(-5) 

2.64(-2) 

1 

1 

0 


1 

1 


2 

0.56 

217304.927 

2.81(-5) 

4.66^2) 

1 

1 

0 


1 

1 


1 

0.33 

217277.680 

2.76(-5) 

3.00(-2) 

1 

1 

0 

<— 

1 

1 


0 

0.11 

217264.639 

2.76(-5) 

1.00(-2) 

1 

1 

1 

<— 

1 

1 


2 

0.10 

217072.801 

5.01(-6) 

8.29(-3) 

1 

1 

1 


1 

2 


2 

0.31 

217074.239 

1.55(-5) 

2.57^2) 

1 

1 

1 

4 - 

1 

1 


0 

0.08 

217032.603 

2.00(-5) 

7.27(-3) 

1 

1 

1 

4 — 

1 

2 


1 

0.10 

217046.988 

8.34(-6) 

9.06(-3) 


Table 2. Integrated intensities (Tdv mK km s main-beam temperature scale) from Gaussian fits to hfs lines N = 0 ■<— N' = 1, 
Fi = 1, observed in the four outflow sources listed, at the offsets {dx", dy") indicated. The error bars are Icr . 


F2^F' 

F <- F' 

u MHz 

L1544(0, 0) 

L183(0, 0) 

L1517B(-10, -20) 

Oph D(0, 0) 

1 -i- 2 

2 <- 3 

108780.201 

73.5 ±3.9 

32.0 ±3.0 

31.3 ±2.9 

25.7 ±6.8 

1 ^ 2 

1 2 

108782.374 

43.2 ±3.8 

12.0 ±3.0 

13.7 ±2.0 

11.9 ±2.5 

1 <- 2 

0^1 

108786.982 

19.3 ±4.0 

10.0 ±2.0 

9.3 ± 3.5 

6.5 ±2.1 

1 <- 2 

1 1 

108793.753 

19.3 ±4.5 

0.0 ±0.0 

7.6 ±2.4 

5.1 ±3.6 

1 ^ 2 

2^2 

108796.400 

18.9 ±3.3 

0.8 ±0.9 

5.2 ± 2.0 

6.0 ±3.1 
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Table 3. Relative intensities of hfs lines = 0 ■<— N' = 1, Fi = 1, observed in the four starless cores listed, at the offsets {dx", dy") 
indicated. The values calculated from the results in Table 1 are also given. 


F2^F' 

F 4 - F' 

V MHz 

L1544(0, 0) 

L183(0, 0) 

L1517B(-10, -20) 

Oph D(0, 0) 

Calc. 

1 2 

2 4-3 

108780.201 

1.00 

1.00 

1.00 

1.00 

1.00 

1 -i- 2 

14-2 

108782.374 

0.59 ±0.06 

0.38 ±0.10 

0.44 ± 0.08 

0.46 ±0.16 

0.54 

1 -i- 2 

0 4- 1 

108786.982 

0.26 ±0.06 

0.31 ±0.07 

0.30 ±0.11 

0.25 ±0.10 

0.25 

1 2 

1 4- 1 

108793.753 

0.26 ±0.06 

0.00 ±0.00 

0.24 ±0.08 

0.20 ±0.15 

0.19 

1 2 

2^2 

108796.400 

0.26 ±0.05 

0.03 ± 0.03 

0.16 ±0.07 

0.23 ±0.14 

0.18 


= ( 11 ) 

where m is the mass of the electron and is the line 
profile function; the spontaneous radiative transition prob¬ 
ability and the absorption oscillator strength are related by 


LJjAji 


mc^ 


^ifij 


( 12 ) 


(Nicholls & Stewart 1962). If the line profile is taken to be 
Gaussian, then 


uo Vtt/ 


“•5 Piv-vof 

exp- 2 -> 


(13) 


where j3 = MF/{2kTj:>), M is the mass of the molecule 
and Td is the corresponding Doppler temperature; vq is the 
frequency at the centre of the line. If Ai/ is defined as the 
(1/e) half-width of the line, then (/(t'o) = l/(-^/7rA!^). 

Let us consider the strongest of the hyperfine 

transitions, at a frequency of 108780.201 MHz. From (12) 
and Table 1, the absorption oscillator strength of this tran¬ 
sition is fij — 1.68 X 10~®. Integrating equ. (11) over the line 
of sight, and assuming that the fractional level population, 
n(i”cN) ’ constant, the optical depth at the line centre is 


t{uo) 


nC^cn) 


Hi -kF 

n(i®CN) me 





where A(^^CN) is the column density of the molecule. 

In L1544 and L183, Hily-Blant et al. (2008, tables 1 and 
2) list column densities V(^®CN) « 10^^ cm“^ at offset (0, 
0); from their figs. C.l and C.2, [S.v « 0.2 km s“^, equivalent 
to Ai/ « 0.07 MHz at the line frequency. Our calculations 
yield ni/n(^®CN) « 10~^. It follows that t{uq) « 10“^. 
For column densities A(^^CN) ^ 10^"^ cm“^, r(i/o) ^ 1, 
and deviations from the LTE line ratios are to be expected. 
Departures of hyperfine level populations from LTE have 
been attributed to line-optical-depth effects in the cases of 
N 2 H'*' (Daniel et al. 2006) and, possibly, NH 2 D (Daniel et 
al. 2013). 

In the following Section 3, we assume the hyper¬ 

fine transitions to be optically thin and calculate their in¬ 
tensities following Section 2.3. 


3 RESULTS AND DISCUSSION 

The relative level populations and line intensities have been 
computed for a range of values of the kinetic temperature. 


T, and the para-H 2 density, n(H 2 ). In Table 1 are listed the 
values of the relative line emission, I [equ. (10)], determined 
separately for transitions within the A = 0 <— A' = 1 and 
A = 1 ■«— A' = 2 rotational manifolds and for T = 10 K and 
n(H 2 ) = 10® cm“®, values typical of pre-protostellar cores; 
the CBR was included in these calculations. 

The values of the relative line emission in Table 1 for 
A = 0 ■4— A' = 1 differ insignificantly from the correspond¬ 
ing results in table A.l of Hily-Blant et al. (2008), calculated 
in LTE - i.e. assuming that the populations of the hyper¬ 
fine states of given A are in proportion to their statistical 
weights, (2F -\- 1). In fact, we found negligible variations of 
the relative line emission within a rotational manifold as 
n(H 2 ) was varied from the low to the high density limit, in¬ 
cluding either the CBR or a black-body radiation field that 
is assumed to be emitted by dust grains in and at the same 
temperature as the gas. This result reflects the fact that 
neither collisional nor radiative transitions discriminate be¬ 
tween the hyperfine states, F, with a given value of A, and 
that each magnetic sub-state, Mf, is equally populated - 
a rare example of LTE conditions prevailing in the ISM. 
Whilst it has long been assumed that LTE applies to the 
relative populations of molecular hyperfine levels (cf. Walm- 
sley et al. 1982), the present study is the first to demon¬ 
strate the validity of this assumption, in the case of ^®CN. 
An analogous conclusion was reached by Daniel et al. (2006) 
regarding N 2 H''', in the optically thin regime. 

Unlike the hyperfine transitions of CN and HCN, which 
have been considered by Faure & Lique (2012), the lines of 
are likely to be optically thin (cf. Section 2.4), given 
that the isotopic abundance ratio « 1:70 in the ISM 

(Milam et al. 2005). It is surprising, therefore, that Hily- 
Blant et al. (2010) concluded that the relative intensities of 
some of the A = 1 —>■ 0 hyperfine transitions of ^®CN, ob¬ 
served in four starless cores, depart significantly from LTE. 
However, a close scrutiny of their fig. 3 reveals that the rel¬ 
ative “LTE” values in this figure are inconsistent with the 
correct values, given by Hily-Blant et al. (2008, table B.2). 

We have re-analyzed our IRAM 30 m telescope obser¬ 
vations of the four starless cores L1544, L183, L1517B, and 
Oph D, and the integrated line intensities - at the offset for 
which the lines are strongest in each source - are reported in 
Table 2; the intensities relative to the strongest transition, 
at 108780.201 MHz, are given in Table 3, together with the 
corresponding theoretical results from Table 1. Our conclu¬ 
sion from this comparison is that the relative populations of 
the hyperfine levels of within a given rotational state, 

A, do not depart significantly from LTE in these sources. 

Although LTE is found to apply to the relative popu¬ 
lations of hyperfine levels of a given rotational state. A, the 
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T(K) 



n(H^) (cm') 

Figure 1. The relative intensity, 1(217467.150 
MHz)/I(108780.201 MHz), of the strongest of the transi¬ 
tions between N = 2 and N = 1 (217467.150 MHz) and between 
N = 1 and N = 0 (108780.201 MHz); the relative intensity is 
plotted as a function of both temperature, T (upper panel), and 
para-H 2 density, n(H 2 ) (lower panel). 


same is not true of the relative populations of levels belong¬ 
ing to different values of N. This situation is reminiscent 
of that which obtains in molecules such as H 2 , where the 
relative populations of the rotational levels within a given 
vibrational manifold can attain LTE at densities which are 
much lower than are required to thermalize the relative pop¬ 
ulations of different vibrational states. In Fig. 1 is plot¬ 
ted the relative intensity of the strongest of the transitions 
(El = 1; E 2 ^ F 2 , F ^ F') between Af = 2 and Af = 1 
(2 •<— 3, 3 •<— 4) and between N = 1 and A'^ = 0 (1 t— 2, 
2 •<— 3). Results are given as functions of both temperature, 
T, and para-H 2 density, n(H 2 ); LTE is approached at densi¬ 
ties n(H 2 )> 10^ cm“^, and n(H 2 )= 10'° cm“° is effectively 
the high density limit, at which a Boltzmann distribution of 
population obtains. As may be seen from Fig. 1, the relative 
line intensity is much more sensitive to T than to n(H 2 ), 
for the relevant ranges of these parameters, which opens the 
possibility of using this line ratio as a thermometer in star¬ 
less cores. 


4 CONCLUDING REMARKS 

We have computed the relative population densities of the 
hyperfine states that comprise the N = 0,1,2 rotational 
manifolds of '°CN, in the range of kinetic temperature 

5 ^ T 20 K and of molecular hydrogen density 1 ^ n(H 2 ) 
^ 10 '° cm“°. The rate coefficients for collisional transfer be¬ 
tween the hyperfine states derive from the recent quantum 
mechanical calculations of Flower & Lique (2015), which 
relate specifically to '°CN-para-H 2 . Spontaneous radiative 
transitions, and those induced by the CBR, were taken into 
account. We find that the hyperfine states, F, with a given 
value of the rotational quantum number, N, are populated 
in proportion to their statistical weights, {2F -f 1), as ex¬ 
pected in LTE, even at densities that are much lower than 
required to thermalize the populations of levels with differ¬ 
ing values of N. A re-analysis of observations of hyperfine 
transitions A = 1 —>■ 0 of '°CN, observed in a number of 
pre-protostellar sources, confirms that the hyperfine states 
are, indeed, statistically populated. 
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